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1.  Introduction 

In  the  past  few  decades  there  has  been  extensive  effort  to  obtain 
exact  solutions  to  the  Schodinger  equation  for  oscillator  systems  with  time- 
dependent  Hamiltonians.  The  path  integral  formalism  of  Feynman  provides 
a  general  approach  to  quantum  systems.  In  this  theory  one  must  obtain  the 
exact  propagator  associated  with  the  classical  action  of  a  given  system. 
Provided  the  exact  propagator  is  obtained,  the  wave  function  of  the  system 
can  be  readily  calculated.  Even  though  the  solution  of  the  Schrddinger 
equation  is  possible  through  this  method,  there  remains  the  problem  of 
second-quantization,  which  is  important  in  connection  with  construction  of 
the  explicitly  t  ime  -  dependent  invariant  (dynamical  invariant)  and  finding 
the  coherent  states  of  the  system. 

Since  Lewis  and  Riesenfeld'  derived  a  simple  relation  between  the 

eigenstates  of  the  dynamical  invariant  and  the  solution  of  the  Schrddinger 

equation,  several  authors  have  employed  the  dynamical  invariant  method  to 

investigate  systems  with  a  time -dependent  Hamiltonian,  including  the  time- 

dependent  harmonic  oscillator.  The  so-called  Ermakov-Lewis  problem  and  its 

(4) 

generalization  have  been  investigated  by  the  following  four  methods:  (i) 

exact  adiabatic  invariants,  (ii)  time-dependent  canonical  transformation, 
(iii)  Noether's  theorem  and  (iv)  Lie  theory  of  extended  groups.  The 
dynamical  invariant  involves  an  auxiliary  function  which  is  related  to  the 
amplitude  of  the  classical  harmonic  oscillator  and  satisfies  a  nonlinear 
second-order  differential  equation  (auxiliary  equation).  Each  particular 
solution  to  the  auxiliary  equation  determines  the  dynamical  invariant. 
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Coherent  states  were  used  by  Glauber^'*  to  discuss  the  photon 
statistics  of  radiation  fields.  After  that,  they  have  been  widely  used  in 
various  fields  of  physics.  Hartley  and  Ray ^ ^  constructed  coherent 
states  for  a  time-dependent  harmonic  oscillator  on  the  basis  of  the  Lewis 
and  Riesenfeld  theory.  Recently,  Yeon,  Um  and  George^^  obtained  the 
exact  coherent  states  for  a  damped  harmonic  oscillator  with  constant 
frequency.  Quesne^^  examined  the  unitary-operator  coherent  states. 

The  ordinary  coherent  states  may  be  defined  in  alternate,  but  essentially 
equivalent,  ways.  For  example,  these  states  are  defined  as  the  eigenstates 
of  a  destruction  operator  and  are  also  obtained  by  applying  a  unitary 
operator,  consisting  of  destruction  and  creation  operators,  to  the  ground 
state  of  the  system.  The  coherent  states  have  several  novel  properties , 

including  the  minimum  uncertainty  product  in  position  and  momentum. 

(12) 

In  previous  work  we  considered  a  molecular  system  adsorbed  on  a 
dielectric  solid  surface,  modeled  as  a  damped  harmonic  ascillator  driven  by 
an  external  electric  field.  The  induced  dipole  moment  of  the  adsorbed 
molecule  obeys  the  equation  of  motion 

x  +  yx  +  w^x  -  f(t)/raQ  ,  (1-1) 

where  ,  y  and  are,  respectively,  the  mass  of  the  adsorbed  molecule, 
modified  damping  constant  and  modified  frequency  due  to  the  presence  of  the 
solid  surface.  The  dots  denote  the  derivative  with  respect  to  time  t,  and 
f(t)  is  an  external  driving  force  given  by 


f(t)  -  g (u',9)  cos (a>' t  -  <^Q) 


(1-2) 
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where  g(o>',0)  is  the  amplitude  of  the  driving  force,  which  depends  on  the 
frequency  cj'  of  the  incident  field  and  on  the  incident  angle  defined  with 
respect  to  the  normal  direction  to  the  solid  surface,  and  4> ^  is  the  phase 
determined  by  incident  and  reflected  fields.  The  Lagrangian  L  and  classical 
Hamiltonian  H  corresponding  to  the  equation  of  motion  (Eq.  (1-1))  are 


L(x,x. t) 


1  2  2 

rm.w  x  +  xf(t) 
t  u  m 


,7t 


H(x , p , t) 


- 2yt  1  22 

6  +  ?VmX 


xf  (t) 


e 


7t 


(1-3) 


(1-4) 


where  p  is  the  classical  linear  momentum.  Hereafter  we  shall  use  the  units 
h  (Planck's  constant)  -  -  1. 

In  Sec.  2. A  we  review  the  results  of  our  previous  work  for  finding  the 
wave  function  of  the  damped  driven  harmonic  oscillator  within  the  path 
integral  f  ormalism  .  ^  ^  The  second-quantization  formalism  is  presented 
in  Sec.  2.B.  The  appropriate  destruction  and  creation  operators  are 
defined,  and  by  using  these  operators  we  obtain  the  wave  functions  for  this 
system.  In  Sec.  2 . C  we  use  the  dynamical  invariant  method,  constructing  the 
exact  invariant  by  solving  the  auxiliary  equation.  The  coherent  states  are 
obtained,  and  their  properties  are  investigated  in  Sec.  3.  In  Sec.  4  we 
briefly  discuss  the  results  and  present  conclusions. 
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2 .  Three  Methods 
A.  Path  Integral  Method 


The  Hamiltonian  operator  H(t)  corresponding  to  the 
Hamiltonian  (Eq.  1-4)  is  obtained  by  making  the  replacement  p  -*■ 


classical 

I  L. 

i  dx* 


H(t) 


i  sL.  -27c  122  .... 

•2  a  2  e  +  2“mx  -  xf(t> 
dx 


(2-1) 


where  x  and  p  satisfy  the  commutation  relation  [x,p]  -  i.  In  order  to  solve 
the  time -dependent  Schrodinger  equation 

i  ^  *(x,t)  -  H(t)  0(x,t)  ,  (2-2) 

( 2  ) 

we  follow  Feynman's  path  integral  method  and  adopt  a  Gaussian  type 
propagator  K, 


K(x, t; 


x0’0) 


2  'Yt 

exp( -a^(t)x  e 


-  a2(t)xe 


yt/2 


a3(t) 


(2-3) 


where  Aq  is  a  normalization  constant  and  we  have  used  the  notation  Xq  -  x(0) 
for  simplicity.  We  assume  the  external  driving  force  to  be  turned  on  at 
time  zero  (t  -  0),  and  for  t  <  0  the  system  is  described  by  the  damped 

/  g  ) 

oscillator  wave  function.  The  propagator  contains  all  information  about 

the  system  in  the  time  interval  [0,t]  and  satisfies  the  wave  equation 


i  “  -  H(t)  K  . 


(2-4) 
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From  Eqs  .  (2-1),  (2-3)  and  (2-4)  we  obtain  three  first-order  differential 
equations  for  the  coefficients  in  the  propagator  (2-3).  The  solutions  of 

those  differential  equations  yield  the  explicit  form  of  the 

(12,13) 

propagator,  " 


K(x , t ;  Xq,0) 


u>  e7t/2  12 


[2?ri  sin(wt)J 


exp 


2  2 
A^(t)x  +  A2(t)x  -  A3(t)xQ 


A4(t^xoJ ’ 

(2-5) 


where 


Ax(t) 


a2(c)  -  LuA(u'  .»)e7t/2[{^  -  £(t)3  cos(w'  sin(w' 


A3(t)  -  — £(t) 


A4(t)  - 


7t 


i  sin(wt)"  {X  •  A(w'-5)  COS  <«'t>V*l)} 


£(t)  -  ^  +  cot(wt) 

2  I 

r  2  il\2 

w  "  \wm  4  J 

A(w'  ,9)  -  g («',«)/ 

-1 


(7«'  )2+^w^-u'2J- 


2n2 


-  tan 


7<j'/(w_2-w'  2)  |  . 
m 


(2-6) 


w  represents  the  reduced  frequency,  which  is  assumed  to  be  real  throughout 
this  paper.  This  implies  that  we  are  concerned  with  the  underdamped  case, 
although  the  results  can  be  taken  over  to  the  overdamped  or  critical  damping 


case . 
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The  wave  functions  of  the  system  are  readily  calculated  using  the 
formula 


0n(x, t) 


K[x, t; 


x0,°]  0n(xo>O), 


(2-7) 


where  ^(x^.O)  is  the  wave  function  for  a  simple  harmonic  oscillator  at  t  -  0, 


VX0'0) 


1  2 

Hn^x0)  exP  t  ‘  2WX0^ 


(2-8) 


and  is  the  usual  Hermite  polynomial.  Substituting  Eqs .  (2-5)  and  (2-8) 
into  Eq.  (2-7)  and  performing  the  integration,  we  obtain 


^n(x,t)  - 


POlL 


2nn!  /tt* 


exp 


|^-B1(t)x' 


+  B2(t)x 


B3(t) 


i (n+2) cot  1(|(t)) 


H 

n 


D(t) 


{x-E(t)|j  , 


(2-9) 


where 


D(t) 


-  |we7t/;i(t)|2 


E(t)  -  A(u',0)  cos(w't  -  -  <ji^) 

B1(t)  -  \  D2(t) 


{1  +  ir2(t)| 
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)jl  +  iC3(t)| 


B2(t)  -  D  (t)E(tHl  +  iC 3 ( t ) 


B3(C)  -  ^D2(t)E2(t)  jl  +  i£(t)J> 


(2-10) 


r1(t)  - 


sin2 (ut)  |l  +  f2(t)| 


r2(c) 


-  <«=) +  •  «>=>} 


0) 


r3(t)  -  f2(t)  -  r1(c)  —  tan(o)'  t  -  4>0  -  <A1) 


Hereafter,  we  shall  use  the  notation  D  -  D(t),  E  -  E(t),  f3  ” 
f^(t)  and  f  -  £(t)  whenever  there  is  no  ambiguity.  We  note  that  there  exist 
useful  relations  between  coefficients  in  the  exponential  terms  in  the  wave 
function. 


Bx  +  B*  -  D2,  B2  +  B*  -  2D2E,  B3  +  B*  -  D2E2 


(2-11) 


which  will  be  used  in  later  calculations. 


B.  Second-quantization  method 

To  implement  the  second-quantization  formalism  for  the  damped  driven 
harmonic  oscillator,  we  introduce  the  time -dependent  annihilation  operator 
a(t)  and  creation  operator  a^(t) 


a(t) 


& 

<j> 

m 


u(t) 


{ 


2BlX 


+  ip 


(2-12) 
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af(c)  - 


H 


u  (tH2B.x 


where 


u(t) 


u0(t) 


e^(t) 


uo(t) 


O) 

m 


2D(t) 


y(t) 


-l 

cot 


(2-13) 


x  and  p  are  the  canonically  conjugate  coordinate  and  momentum  such  that 
[x,p]  -  i,  and  B^ ,  B^  and  D  are  all  given  in  Eq.(2-10).  It  is  obvious 
through  Eq.  (2-11)  that  the  non-Hermitian  operators  a(t)  and  a^(t)  satisfy 
the  commutation  relation 


[a(t) ,  af ( t) ]  -  1 


(2-14) 


Thus  we  set 


a(t)|n,t>  -  ,/n|n-l,t> 


a^(t)|n,t>  -  7n+l | n+1 , t> 


a(  t)  |  0 ,  t>  -  0 


(n  -  0,1,2, 3.  .  .) 


(2-15) 
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The  state  vector  |n,t>  implies  that  the  quantum  state  depends  on  the  quantum 
number  n  and  time  t.  Even  though  the  definitions  of  a(t)  and  a'(t)  in  Eq. 
(2-12)  are  correct,  it  is  more  convenient  to  introduce  a  new  function  for 
later  calculations. 


f?(t)  -  —  u(t)  B.  e‘7C 
to  1 

m 


(2-16) 


Then,  the  time-dependent  operators  a(t)  and  a^(t)  are  rewritten  as 


a(t)  -  -  „  2  i  e7C  r,(t)(x-E)  -  e'7t  jp  + 

-  m  ^  - 


(2-17) 


a!(t)  -  -/2i  e7t  e‘7t  |p  +  f4(t)l  -  i,*(t)(x-E) 

-  m  '  ' 


where 


r4(t)  -  D  E  (f2-f3) 


(2-18) 


From  now  on  we  use  the  notation  a  -  a(t),  a^  -  a^(t),  u  -  u(t),  etc.  x  and 
p  are  expressed  in  terms  of  a  and  a^  as 


J1  *f, 

u  a  +  ua 1 


x-E+  (ua+ua) 
u> 

m 


_  /7T  J  t  *  a  a  t 

p  -  +  y2  e  [rj  a  +  } 


(2-19) 
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Substitution  of  Eq .  (2-20)  into  (2-1)  gives  the  second- quantized  expression 
of  the  Hamiltonian  operator  H(t), 


H(t)  -  e 


7 t 


/I  I  2  1  t  2  .  ,  a  t  A  f  2  2  .  a  f  A  j 

( I U I  +  j  17 1  )(2a  a  +  1)  +  ju  -t-rj  )a'a 


*  (^“»e  ■  f<»i  ■  ”^yyc 

m 


)at  +  h-c-}] 


+  G(t) 


G(t)  -  e 


7t 


1  22  ^  12  -27t  .... 

E  +  e  -  Ef ( t) 

l  m  1  4 


(2-20) 


where  h.c.  implies  the  Hermitian  conjugate. 

The  ground-state  wave  function  in  the  coordinate  representation, 
^(x.t)  -  <x|0,t>,  is  readily  calculated  using  the  definition  of  a  in  Eq . 
(2-12)  and  a|0,t>  -  0.  The  normalized  form  is 


Vx,c)  “ 


'  D  | 

2 

IM 

exp 

-  B^x 


BjX 


-  B. 


-  \  y(t) 


(2-21) 


which  corresponds  to  Eq.  (2-9)  when  n  -  0.  The  n-th  state  wave  function 
^(x,t)  is  obtained  through 


*n(x,t)  “  y^f(at)n  VXlt)  • 


With  a^  as  defined  in  Eq .  (2-12),  we  readily  obtain 


(2-22) 
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fe-iyUKn 

1  v'2D  J 


exp 


r  *  2 


exp 


*  2 
■B-jX 


* 

+  B?x 


(2-23) 


Substituting  Eqs  .  (2-21)  and  (2-23)  into  (2-22)  and  performing  some 

calculations  using  Eq .  (2-11),  we  obtain 


tf>n(x,t)  - 


,  \n 

■1)  exp 


-B^x^+B^x-B^- i(n+^)y(t) 


2  a  2 
z  -z 

’  (a2>  e 


(2-24) 


where  z  -  D(x-E) ,  and  the  Hermite  polynomial  function  is  expressed  as 


Hn(z) 


(-l)n  e: 


(2-25) 


Then  the  wave  function  in  Eq.  (2.24)  is  exactly  the  same  as  Eq .  (2-9).  It 

is  easy  to  show  that  (x,t)  has  the  orthonormal  property, 

<V>(x,t)|i/>(x,t,)>-<m,t|n,t>-5  (2-26) 

m  n  m,  n 

C.  Time -dependent  operator  method 

For  the  system  characterized  by  the  time-dependent  Hamiltonian  H(t) 

A 

[Eq.  (2-1)],  we  assume  that  there  exists  a  Hermitian  operator  I(x,p,t)  which 
is  explicitly  time -dependent  and  satisfies  the  invariant  condition 


d_ 

dt 


I 


[ I ,H]  -  0 


(2-27) 


By  operating  on  the  left-hand  side  of  Eq.  (2-27)  with  the  time  -  dependent 
Schrodinger  state  vector  |>  which  satisfies  the  wave  equation 
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1 


H(t)|>s 


(2-28) 


we  obtain 


L- 

3t 


(I 


>s)  "  H ( t ) < I | >s ) 


(2-29) 


This  means  that  the  operation  of  the  invariant  on  the  Schrodinger  state 

(14) 

vector  yields  another  solution  of  the  Schrodinger  wave  equation.  Now, 

(2) 

we  assume  that  the  time -dependent  invariant  operator  has  the  form 


I(x,p,t)  -  - 


61(t)x 


<$2 (t)  (xp+px) 


<53(t)p 


64(t)x 


+  5^(t)p  +  6g(t) 


(2-30) 


where  the  S ^ ' s  are  all  real  functions  of  time  and  do  not  involve  time- 
differential  terms.  Applying  Eq.  (2-30)  to  (2-27),  we  have  six  differential 
equations  for  the  5^'s, 


L  -  2 5_  w2  e7t  -  0  (2-30-1) 

i  z  ra 

5-  +  6.  e‘7t  -  5,  w2  e7t  -  0  (2-31-2) 

l  i  i  m 

53  +  2 S2  e'7t  -  0  (2-31-3) 

5.  +  (25,f  -  5 cw2)  e7t  -  0  (2-31-4) 

4  Zorn 

L  +  6.  e‘7t  +  25, f  e7t  -  0  (2-31-5) 

D  4  3 
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8,  +  {  f  eyt  -  0 

O  J 


(2-31-6) 


where  the  dots  represent  partial  derivative  with  respect  to  time. 
To  solve  the  above  equations,  we  write 


«3(t) 


P2(t) 


-7t 


(2-32) 


where  p(t)  is  a  real  function  of  time  which  will  be  determined  later. 
Substituting  Eq.  (2-32)  into  Eqs .  (2-31-3),  (2-32-2)  and  (2-31-1),  we  get 


5x(t) 


7  P 


2YI 

u>  ► 

mJJ 


,7t 


«2<t) 


(2-33) 


where  p(t)  obeys  the  nonlinear  differential  equation 


p  + 


2 

w  p 


-3 


(2-34) 


and  w  is  the  reduced  frequency  defined  by  Eq .  (2.6).  Any  particular 

A 

solution  of  Eq.  (2-34)  can  be  used  to  construct  the  invariant  I(x,p,t).  It 
is  straightforward  to  show  that 


Pit) 


/2 


/D(t) 


(2-35) 


is  a  particular  solution  of  Eq .  (2-34),  where  D(t)  and  f^(t)  are  defined  in 
Eqs.  (2-10).  Inserting  this  solution  into  Eqs.  (2-31),  we  obtain 
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5x(t)  -  d2[i  +  r2] 


s2(t)  -  r2  -  2<j  '7ri  ' 


(2-36) 


Substitution  of  Eq  .  (2-36)  into  Eqs .  (2-31-4)  and  (2-31-5)  yields  a 

nonhomo geneous  differential  equation  with  constant  coefficients, 


6*5  +  765  +  4*5 


2rtf  +  {7^  +  3rjf  /» 


(2-37) 


the  solution  of  which  is 

fi5(t)  -  -2E  C3  , 


(2-38) 


where  E(t)  and  c )  have  been  defined  in  Eqs.  (2-10).  Again,  substituting 
Eq.  (2-38)  into  (2-31-4),  (2-31-5)  and  (2-31-6),  we  obtain 


54<t)  -  -2D2E  (1  +  ;2$-3) 
2  2  2 

*6(t)  -  dV  (i  +  f‘)  . 


(2-39) 


We  have  determined  all  coefficients  in  Eq.  (2-30).  Thus  the  invariant 

A 

operator  I  can  be  explicitly  written  as 


I(x,p,t) 


If  2  2  (  p  11 

2  D^(x-Er  +  |D(t)f2x  +  -  -  DEf3j 


(2-40) 


To  obtain  the  operator  form  of  I ,  we  introduce  the  two  operators 
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aQ(t)  -  ^  uQ(t)  (2B1x  +  ip  -  B2) 
m 

(2-41) 

*0(t)  "  ^  u0(t)  (2B1X  '  ip  •  B2!  ’ 
m 


where  Ug(t)  is  defined  in  Eq .  (2-13).  We  note  that  ag(t)  -  a(t)e  and 
aj(t)  -  a^tje^^^  .  Equation  (2-40)  is  expressed  in  terms  of  ag(t)  and 
aj(t)  as 

I(t)  -  aJ(t)aQ(t)  +  i  .  (2-42) 


Obviously,  a^(t)  and  a^(t)  satisfy  the  same  commutation  rule  as  a(t)  and 
a*(t),  i.e., 


[aQ(t),  aj(t)]  -  1 


(2-43) 


Therefore,  we  can  write 


aQ(t)|n,t>I  -  /n| n- 1 ,  t>z 

a^(t)jn,t>1  -  yn+1 | n+1 , t>1  (2-44) 

a0(t)|n,t>z  -  0  (n  -  0,1,2,3,...) 

A 

Here  we  assume  that  the  eigenfunctions  ^  (x,t)  -  <x|n,t>z  of  I  form  a 
complete  orthonormal  set  corresponding  to  the  eigenvalues  X Then  we  have 
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A 

16  (x.t)  -  X  6  (x.t) 
n  n  n 

A  -  n  +  i  .  (2-45) 

n  L 

A 

Even  though  the  eigenfunctions  of  I  are  t ime- dependent ,  the  eigenvalues  are 

(3) 

time -  independent . 

In  order  to  obtain  6  (x.t),  we  first  obtain  the  normalized  form  of 

n 

.^(x.t)  using  the  explicit  expression  of  £g(t)  in  Eq.  (2-44), 

1 

0Q(x,t)  ”  2  exP  +  B2x  -  B3  ,  (2-46) 

which  is  the  same  as  except  for  the  y(t)  term  in  Eq.  (2-21).  The 

eigenstate  <^(x,t)  can  be  obtained  by  applying  a^t)  n  times  successively  to 
^(x.t).  Through  the  same  procedure  for  obtaining  Eq.  (2-24),  we  get 

*.(>•'>  "if  <4,n*o(It’c) 

1 

(2-47) 

Now  the  solution,  V’n(x,t),  of  the  Schrodinger  wave  equation  can  be 
obtained  fron/2^ 

ia(t) 

^  (x,t)  -  e  d  (x,t)  ,  (2-48) 

n  n 

where  the  phase  function  a^t)  is  the  solution  of  the  equation 
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(t)  o 

— T -  =  <4>  |  i  -  H  U  > 

dt  n'  ot  1  n 


(2.49) 


The  diagonal  matrix  element  of  —  with  respect  to  the  eigenvector  |n,t>T  of 

Ot  -L 

I  is  obtained  by  taking  the  partial  time  derivative  of  a^(t)  i-n  Eq-  (2-44) 
and  operating  taking  a  scalar  product: 


1<n,tj-|n,t>1 


r<n- 1 , t 


Sal 


-in-l.t^  +  ^  j^.tl 


—  Jh-I, 


a 

[<0  ’  1 1  a"t  I  0  ’  t>i  +  ^  I<n.  t|  —  In-i.t^ 


(2-50) 


In  obtaining^  the  second  line  in  Eq .  (2.50),  we  have  used  a  r 

al7 


ecurrence 


relation.  — “  is  easily  calculated  from  Eq.  (2-41)  as 


d&0  2i  yt  /  *  •  F  2i  ytf  •*  ,  •  xllgt 

~t  -  “  e  t?7o(7Uo-uo)  +  VoJao  +  7  +  w  e  VVo  +  ,Jo(7Uo-uo)|Jao 

ra  v  1  L  m  v 


U0'4  +  “m  " 


;  * *1 


(2-51) 


where  q g(t)  is  a  function  which  is  obtained  from  Eq.  (2-16)  by  making  the 
replacement  u(t)  -*  UQ(t).  Thus,  Eq.  (2-50)  becomes 

I<n  ’  1 1  ft  t n .  t>x  -  ^O.tlf^lO.Oj  +  n[7  +  e7t{u0^  +  V7VV}]  . 

L  m 

(2-52) 


^<0,  t|^|0,  t>j  gives  the  zero-point  contribution  to  Eq.  (2-52),  which  is 
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-  ^[7  +  ^  {uQ^o  +  V^VV}] 

L  m  v 


(2-53) 


Using  Che  definition  of  uq(c)  and  r/^t)  and  Eqs .  (2-10),  we  obtain 


Vo 


V^W 


Tw  e 
2  m 


-7t 


2  r2 


(2-54) 


Substitution  of  Eqs.  (2-53)  and  (2-54)  into  (2-52)  yields  the  diagonal 
matrix  element  of  i 

o  t 


i<n, c I i  ^!n-c>x 


(2-55) 


A  A 

Now  we  can  express  the  Hamiltonian  H  in  terms  of  ag(c)  and  ag(t)  by 
making  the  replacements  u(t)  -*  Ug(t)  ,  rj( t)  -*  nQ(t)  ,  a(t)  -*  aQ(t)  and  a^(t)  ■+ 
aj(t)  in  Eq.  (2-20).  By  this  manner  we  obtain  the  diagonal  matrix  element 

A 

of  H  as 


j<n, t|H|n, t>x 


2(n  +  2)  (uo  +  l’7ol2)e7t  +  G(t) 


(2-56) 


Here,  we  note  that  the  G(t)  term  in  Eqs.  (2-20)  and  (2-56)  does  not  have  any 

effect  on  the  dynamics  of  the  system.  Therefore  it  is  always  possible  to 

remove  this  term.  In  later  calculations,  we  neglect  this  term. 

u 

Since  p  -  (f^/w)  i-s  solution  of  the  auxiliary  equation  (2-34), 
this  equation  can  be  written  in  the  equivalent  form 


w  + 


2 

X_ 

8w 


2WC1 


(2-57) 


In  obtaining  the  above  equation  we  have  used  the  identity 


2 

•*  2  2 

ri  “  +  2  '  ri  '  (2-58) 

which  is  easily  verified  by  using  Eqs .  (2-10).  With  use  of  Eq.  (2-57)  we 
obtain 


w  2 


(2-59) 


Using  these  results,  we  finally  obtain 


d^t) 

dt 


(2-60) 


To  integrate  the  above  equation,  we  use 


y(t)  -  i  , 


(2-61) 


which  is  verified  by  using  the  definitions  of  C^(t)  and  y(t).  Thus  we 
obtain 


<*n(t)  -  (n  +  cot  ^[£(t)]  +  constant 


(2-62) 


where  the  constant  gives  a  constant  phase  in  the  wave  function  (x,t),  and 


n 


can  thus  be  neglected.  Substitution  of  Eq.  (6-62)  into  (2-48)  yields  the 
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Schrodinger  wave  function,  which  is  exactly  the  same  as  Eqs .  (2-9)  and  (2- 
24).  Thus,  we  have  shown  that  the  three  methods  described  in  this  section 
yield  the  same  wave  function. 

3.  Coherent  states 

In  this  section  we  obtain  the  coherent- state  wave  function  by  means  of 
three  equivalent  definitions  and  show  that  the  coherent  states  correspond  to 
minimum  uncertainty  states.  We  define  a  coherent  state  | a , t>  as  the 
eigenstate  of  a(t)  corresponding  to  the  eigenvalue  a, 

a ( t ) | a , t>  -  a | a , t>  ,  (3-1) 

where  a  is  a  complex  number  and  a(t)  has  been  defined  in  Eq .  (2-12)  or  (2- 

17).  Equation  (3-1)  can  thus  be  rewritten  in  the  coordinate  representation 
as 

^  u(t)/2B^(t)x  -  B2(t)  +  j-  <x | a ,  t>  -  a  <x|a,t>  (3-2) 

m  ^  * 

^a(x, t)  S  <x|a, t> 

-  N(t)  expj^-B^x2  +  |b2  +  y=^j-  xj  ,  (3-3) 

where  N(t)  is  a  normalization  factor.  Through  the  normalization  procedure 


we  obtain 
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where  we  have  used  the  relation  u/u  =  e2i-y(c)  =  exp[cot  r  ( t )  ]  .  The  wave 

function  0  (x,t)  can  be  obtained  through  the  second  definition  of  the 
a 

coherent  states  such  that 

0  (x,t)  =  exp  [  aa  1  ( t )  -aa(t)Jv>rt(x,t)  (3-5) 

a  '  u 

where  ^(x.t)  is  the  ground-state  wave  function  given  by  Eq .  (2-21).  Thus 

1  'A* 

exp [ aa'  -  a  a]  plays  the  role  of  a  displacement  operator.  If  we  use  the 
identity 


At  i  1 2  . „ 

aa  -q  a  aa1  -a  a  -a  /2 
e  —  e  e  e 


(3-6) 


equation  (3-5)  can  be  rewritten  as 


oo 

,  .  ,  -lal^/2  \  a^  k  r  _\ 

0a(x,t)  -  e  /  n!  (a  ^  V^Cx.t) 


n-0 


-  e 


-  «r/2 


I  *  V- 


t)  . 


n-0 


(3-7) 
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where  we  have  used  Eqs .  (2-15)  and  (2-22).  Eq .  (3-7)  can  also  be  used  for 
Che  definition  of  the  coherent  states.  Since  i/>^(x,t)  is  given  by  Eq .  (2-9) 
or  (2-24),  Eq .  (3-7)  becomes 


V  (x, t)  = 

r  2 

fD  u  1  4 

I  exp 

L  ™  J 

r  2 

- B, x  +  B,x  -  B„ 

L  15 

1|  |2l 

-  2  i a  ' 

z 

,  2 

1  fa  u  ln/2 

n!  [  2u  J 

n»0 

X 

Hn[D(x-E)] 

. 

(3-8) 

The  summation  in  Eq .  : 3-8)  can  be  easily  carried  out  bv  use  of  the  expansion 

of  the  Hermite  polynomial  H  : 


n! 


2  ~ 
ro  u 


2u 


Hn[D(x-E) 


-  exp 


n-0 


ttD(x-E)  J 2u  /U  - 


2  *■ 
a  u 


2u 


(3-9) 


Substitution  of  Eq.(3-9)  into  (3-8)  yields  Eq .  (3-4).  Thus  we  have  shown 
that  three  definitions  of  coherent  states  [Eqs.  (3-1),  (3-5)  and  (3-7)]  are 
equivalent.  Here  we  note  that  when  there  is  no  driving  force  (i.e.  ,  f  -  0, 
E  -  0) ,  the  coherent  state  wave  function  [Eq.  (3-4)]  is  reduced  to  that  of 
the  damped  (not  driven)  harmonic  oscillator , 


0  (x, t) 

a 


r  2  *- 

Du 

U 

n\x 

exp 

-B  x2  +  aD  7— 
1  u 


2 

a  u 
2u 


(3-10) 


We  now  show  that  the  coherent  state  vectors  form  a  nonorthogonal 


complete  set: 
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easily  obtain 
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\0.t>  -  ;  |  d2a  | a , t>  exp  -^(|o|2  +  \fi\2)  +  a* 0 


(3-13) 


This  means  that  Che  coherent  states  are  not  linearly  independent  of  one 
another.  Even  though  these  states  form  a  complete  set,  there  are  more 
states  than  are  necessary  for  expanding  any  given  state  in  terms  of  the 
coherent  states.  In  this  sense,  these  states  are  said  to  be  " overcomplete" . 

The  coherent  states  correspond  to  the  minimum  uncertainty  states,  and 
thus  we  can  show  easily  the  following  relation: 


((Ax)(AP)]q  =-  £(1  +  f2)4 


(3.14) 


Equation  (3.14)  is  the  minimum  uncertainty  corresponding  to  the  ground 
state . 

We  now  show  that  the  coherent  states,  which  are  eigenstates  of  the 
destruction  operator  defined  by  Eq.  (2-12)  [or  Eq.  (2-17)],  give  the  exact 
classical  motion  for  the  damped  driven  oscillator  system.  The  calculation 
of  the  position  of  the  wave  packet  yields 


<a , t | x | a , t> 


-7t/2 


*  iwt  ^  *  *  -iwt 

F(a,a  )  e  +  F  (a, a  )  e 


g(w’  ,  9  )cos(w' t-^-^) 

+  /,  2  ,2,2  .  ,  .  2  ~ 

y(w -w  )  +  (yw') 

m 


F(a,a  ) 


(3-15) 


where  we  have  used  the  explicit  expression  for  u(t)  from  Eq .  (2-13), 
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.  ^m  -yt/2  ( 

u(t)  -  iyS  «  { 


y  .  .  .  iwtl 

2w  sin(ut)  +  e  | 


which  is  readily  obtained  by  use  of  the  identity 


(3-16) 


exp 


±in 


cot'L(£(t)) 


£ ( t)±i 

,£(c)+i. 

(3-17) 


where  £(t)  -  ~~~  +  cot(cjt)  [see  Eq.  (2-6)].  On  the  other  hand,  the  classical 
solution  of  the  equation  of  motion  for  the  damped  driven  harmonic  oscillator 

[Eq.  (1-1)]  is  given  by 


xx(t)  -  e 


-7t/2 


iwt 

e  +  A^  e 


•iwtj 


g(o' , ff)cos(w' t -<i>0-4>l) 

7(W2-W'2)2  +  (7c' 5 


(3-18) 


4.  Discussion  and  conclusions 

We  have  considered  a  molecule  adsorbed  on  a  dielectric  solid  surface, 
modeled  as  a  damped  harmonic  oscillator  forced  by  a  time -dependent  electric 
field.  The  exact  wave  functions  of  this  system  have  been  calculated  by  - 
three  different  methods.  These  methods  yield  exactly  the  same  results. 
When  there  is  no  external  perturbing  force,  the  wave  functions  naturally 
reduce  to  those  of  a  damped  harmonic  oscillator.  Furthermore,  when  there 
exists  no  damping,  the  wave  functions  become  the  time  -  dependent  wave 
functions  for  a  simple  harmonic  oscillator.  In  this  sense,  we  can  say  that 
we  have  obtained  the  general  wave  function  for  a  harmonic  system. 

The  Hamiltonian  in  Eq.  (2-20)  has  off-diagonal  terms  since  the 
representation  defined  by  Eqs.  (2-15)  does  not  diagonalize  this  operator. 
When  there  is  no  external  driving  force,  Eq.  (2-20)  becomes 
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H  -  e 


7t 


/I  I  2  |  |  2.  A  .  .  .  2  2  .  A  t  A  t 

( I  u  |  +1*71  )(2a'a  +  1)  +  (u  +rj  )a  a1  + 


,  *2  *2, 

(u  +r)  ) 


A  A  | 

aa 


(4-D 


The  above  equation  has  the  most  general  form  for  a  Hamiltonian  which 

(16) 

preserves  the  coherence  of  an  arbitrary  initial  SU(1,1)  coherent  state. 

Thus  this  Hamiltonian  can  be  used  to  describe  the  dynamics  of  the  SU(1,1) 
coherent  states  of  a  damped  oscillator.  When  there  is  no  damping  the 
Hamiltonian  in  Eq.  (4-1)  becomes  that  of  a  simple  harmonic  oscillator.  In 

i  iwoc  i  iwoc 

this  case,  we  easily  obtain  u(t)  -  e  ’  “  2^0  e  '  ^*en  we 

A  *  I 

have  H  —  w_(a,a+~).  Therefore,  we  can  say  that  Eq .  (2-20)  is  the  general 

form  for  the  Hamiltonian  of  a  (damped,  driven)  harmonic  oscillator  system. 

To  obtain  the  explicitly  time-dependent  (Hermitian)  invariant 
operator,  we  have  assumed  the  form  of  the  invariant  to  be  expressed  by  Eq. 
(2-30).  The  same  type  of  invariant  was  used  by  Khandekar  and 
Lawande^^  for  a  damped  driven  oscillator  with  a  time -dependent  frequency. 
When  the  frequency  of  the  system  depends  on  time,  the  auxiliary  equation 
necessary  for  constructing  the  invariant  should  be  solved  numerically.  For 
the  constant  frequency  case  like  ours,  the  auxiliary  equation  (2-34)  is 
easily  solved.  The  particular  solution  to  Eq .  (2-34)  is  given  by  Eq.  (2- 
35).  With  this  solution  we  obtain  the  explicit  form  of  the  dynamical 
invariant,  Eq.  (2-40).  When  there  is  no  driving  force,  Eq.  (2-40)  can  be 
written  as 


C  i|",,X.2  7 1  /  1  k  7t]2  -7t 

1  -  j  <p)  e  +  1PP  -  *(P  ’  2p>e  J  e 


(4-2) 


This  equation  agrees  with  the  result  of  Pedrosa^^  and  is  somewhat 

/  IQ  \ 

different  in  form  from  that  of  Korsch.  Here,  we  point  out  that  for  a 

given  equation  of  motion  there  can  be  various  kinds  of  auxiliary  equation*. 


2 


Hence,  a  given  Hamiltonian  can  have  many  different  dynamical  invariants.  Of 
course,  the  same  physical  results  must  be  obtained  regardless  of  which 
dynamical  invariant  is  chosen.  For  the  case  of  no  damping,  Eq .  (4-2) 

(19) 

becomes  the  so-called  Ermakov-Lewis  invariant  for  a  harmonic  oscillator, 


:  ix2. 

1  -  7  ~  +  (pp 

^UpJ 


-xp)2}  . 


(4-3) 


The  operators  ap(t)  and  ag(t)  in  Eq.  (2-41),  which  were  used  to  diagonalize 

A  A  * 
the  invariant  I  [Eq.  (2-42)],  are  the  same  as  a(t)  and  a' (t)  in  Eq.  (2-12) 

except  for  the  phase  factor.  Therefore,  the  phase  function  a^(t)  in  Eq . 

(2-49)  is  readily  obtained  with  use  of  Eq.  (2-20). 

The  coherent  states  defined  as  eigenstates  of  the  operator  a(t)  form  a 

nonorthogonal  (over)  complete  set  and  correspond  to  the  minimum  uncertainty 

states.  The  coherent  states  can  be  defined  in  equivalent  ways  as  described 

in  Sec.  3.  The  product  of  uncertainty  in  position  and  momentum  has  a 

periodic  characteristic  with  period  x/w.  The  uncertainty  product  in  the 

damping  oscillator  system  is  not  affected  by  the  external  perturbing  force. 

For  the  7  -  0  case,  the  uncertainty  product  has  a  minimum  value  of  The 

coherent  states  obtained  satisfy  the  Schrddinger  equation  and  give  the  exact 

classical  motion  of  the  damped  driven  oscillator  system. 
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